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Abstract. We consider a spatially-heterogeneous generalization of a well- 
established model for the dynamics of the Human Immunodeficiency Virus-type 
1 (HIV) within a susceptible host. The model consists of a nonlinear system of 
three coupled reaction-diffusion equations with parameters that may vary spa¬ 
tially. Upon formulating the model, we prove that it preserves the positivity of 
initial data and construct global-in-time solutions that are both bounded and 
smooth. Finally, additional results concerning the local and global asymptotic 
behavior of these solutions are also provided. 

1. Introduction. Over the past few decades, considerable effort has been devoted 
to modeling the in-host dynamics of viral infection, and in particular, HIV in¬ 
fection within humans. These models have made significant contributions to the 
understanding of HIV pathogenesis in vivo and the mechanisms through which the 
infection may be mitigated. The typical response to viral infection within a host 
is the activation of the immune system, driving the level of virions down. If the 
immune response is sufficiently potent then the disease can be completely eradi¬ 
cated from the body, but often this does not occur. Instead, over a time period 
that can vary from weeks to months, an eventual balance of viral replication and 
clearance of the virus by the immune system occurs, leading to a state known as 
chronic infection. These equilibrium outcomes - viral clearance versus the devel¬ 
opment of a chronic infection - are suggestive of simple dynamics, but this ignores 
spatial variations including possible localized regions of infection, which alter the 
behavior in both transient and long time asymptotics [11]. 

Although some mathematical modeling of acute HIV infection has been per¬ 
formed [4, 18, 29], many current models focus on the time course of infection during 
the chronic stage, such as the response to antiretroviral therapy and HIV levels after 
the viral setpoint has been established [1, 27]. Moreover, because of the intricate 
complexity and the enormous number of cells involved, mathematical descriptions 
have generally been limited to nonlinear systems of a few coupled ordinary differ¬ 
ential equations describing the average behavior throughout the whole body under 
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the assumption that the environment is well-mixed or spatially-homogeneous [9, 25]. 
Unfortunately, such an assumption is not valid during the earliest stages of infection 
or at sites of viral entry, for instance, the vaginal or rectal epithelium during sexual 
transmission of HIV. 

Indeed, viral propagation is a fundamentally local process. Focal bursts of viri¬ 
ons have been demonstrated near infected T-cells in vivo [11, 17, 28]. Correlations 
between spatial location and viral genetics have been observed and modeled within 
splenic tissue [3, 7]. Localization is particularly relevant in the case of HIV infec¬ 
tion, for a number of reasons. First, the virus is inherently unstable, possibly due to 
shedding of the glycoprotein gpl20, which mediates the virus-cell binding process 
[15]. Spatially-heterogeneous outcomes may arise from underlying heterogeneity 
within tissue architecture, but it is also possible to observe such non-uniformities 
emerge spontaneously from the infection dynamics. Finally, an important aspect of 
the in-host spread of HIV is that the vast majority of infection occurs only in lym¬ 
phoid tissues, where target cells are densely packed. This particular environment 
allows the virus to maximize the efficiency of diffusive transfer from infected pro¬ 
ducers to target cells. Thus, the local propagation of HIV within lymphoid tissue is 
fundamentally different in comparison to the hematogenous spread of the virus to 
other distant tissues [10, 31]. For all of these reasons, advances in modeling spatial 
aspects of the development and spread of viral infection within a host are crucial 
to furthering our understanding of viral pathogenesis and treatment. 

Though this paper constitutes the first mathematical study of a model using a 
system of partial differential equations to study the spatial dependence of infection 
dynamics, others have used different mathematical mechanisms to model the incor¬ 
poration of spatial fluctuations and correlations. Funk et al [8] posited a discrete 
lattice model to study the spatial dynamics of virions and T-cells, but provided 
only computational results and did not allow for motion or infection away from 
a fixed grid. Strain et al [31] studied the competition between viral lability and 
diffusion using a spatial cellular automaton model with a fixed T-cell background. 
Finally, Brauner et al [2] studied a system of coupled ODEs and a single PDE that 
allowed for the diffusion of virions within a two-dimensional medium under the as¬ 
sumption that T-cells remained motionless. However, their study did not focus on 
the well-posedness or global dynamical behavior of the model, or allow for cellular 
diffusion. 

In order to investigate the impact of spatial dynamics in a simple mathematical 
model of HIV infection we extend the standard lumped or three-component model 
of in-host viral dynamics [1, 13, 19, 20, 26] to include spatially random diffusion 
and a spatially-dependent T-cell supply rate. Upon describing the model, we estab¬ 
lish basic results concerning the well-posedness of smooth solutions and then focus 
mainly on their longtime asymptotic behavior. The overall goals are to establish a 
mathematical framework under which the spatial propagation of infection may be 
studied and to elucidate the contribution of the parameter space on the dynamical 
properties of solutions. 

This paper proceeds as follows. In the next section, we introduce the spatial 
model of in-host viral dynamics and state our main results concerning existence, 
uniqueness, regularity, and longtime behavior of solutions to this system of semi- 
linear PDEs. The proofs of these theorems are contained within Section 3, and 
finally proofs of associated lemmas are included within Section 4. Throughout the 
paper C will denote a positive constant that may change from line to line. When 
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necessary, we will denote the dependence of this constant on other parameters using 
subscripts, for instance Cf ■ 

2. Spatial Model of Viral Dynamics and Main Results. 

2.1. Derivation of the model. We begin with the standard, spatially-homogeneous 
model for in-host virus dynamics, given by 


dT ^ 

— = A - htT - kTV, 
dt ^ 


( 1 ) 



dt 


Here, the unknown quantities are the populations of (i) uninfected target cells, 
denoted by T(t) (in the study of HIV, these are CD4 T-cells); (ii) infected cells, 
denoted by /(t); and (hi) free virions V{t). It is assumed within the model that 
target cells are supplied at a constant rate A and removed either through infection 
via contact with virions at a rate of k per virion, or through natural cell death with 
per capita rate pT- Similarly, fij represents the per capita rate at which infected 
cells are destroyed, either through natural cell death, interaction with the body’s 
immune response, or via lysis due to new virions bursting from the cell membrane, 
while fly represents the per capita rate at which virions are cleared from the body. 
Additionally, N represents the “burst rate” of the virus, or the average number of 
new virions produced over the entire lifespan of an infected cell. 

In order to introduce spatial variations within each population, we first let H C 
R" be a given bounded domain with dfl smooth. As we desire a reformulation that 
preserves the general biological mechanisms of infection dynamics, we consider the 
system 


dtT - DtAT = X{x) - ^J.TT - kTV; 

dtl - DiAI = kTV - fill, 
dtV - DyAV = Nflil - flyV 


(3CM) 


for a; € H, f > 0, which describes a spatially-heterogeneous three-component model 
in which all populations may diffuse throughout a fixed region of the body and the 
supply rate of new susceptible cells is spatially-dependent, rather than constant. 
Hence, T, /, and V now represent concentrations of healthy cells, infected cells, and 
virions, respectively. The diffusive nature of cell transport is, by now, well-known 
both in-vivo and in-vitro, and while there may be some minor evidence for this 
motion to occur via fractional or anomalous diffusion [12], we assume a standard 
diffusive mechanism via Brownian motion as suggested by [16]. As spatial effects 
are present, the domain H might be considered as an entry point of viral infection, 
such as the vaginal or rectal epithelium during sexual contact. Additionally, one 
does not expect the introduction of new T-cells to occur uniformly throughout the 
domain, and thus, we allow for spatial variations in the function A(a;) rather than 
taking this parameter to be constant as in (I). 

With the spatial domain fixed, the coupled, nonlinear system of PDEs is aug¬ 
mented by the initial conditions 


TiO,x) = To{x), I{0,x)=Io{x), ViO,x) = Vo{x). 
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Throughout, we make some standard biological assumptions on the initial concen¬ 
trations, namely 

Tq{x), Io{x),Vq{x) are continuous, positive, and bounded functions on Q (Al) 

Additionally, we consider the separate cases of homogeneous Dirichlet or Neumann 
boundary conditions, so that (3CM) is supplemented by either of the two conditions 


w{t, x) = 0 on [0, oo) X 

O 

x) = 0 on [0, oo) X dfi 
an 


(BCl) 


(BC2) 


for w replaced by T, I, and V, respectively. 

The biological parameters k,N,iJLT,iJii, /iy and diffusion coefficients Dt, Dj, Dy 
are all positive constants. For simplicity, we assume throughout that A € C°° (fl) 
with A(x) > 0 for X € n and A ^ 0. 

With this framework in place, we make a few remarks concerning the generality of 
the model. First, this model can be used to study a wide-range of viral infections, 
including HIV, Hepatitis B (HBV), Hepatitis C (HCV), Human T-cell Leukemia 
Virus (HTLV), and Human Cytomegalo Virus (CMV) [1, 20], but we mainly focus 
on the case of HIV as in-host modeling has significantly furthered the understanding 
of this particular disease in recent years. Additionally, though we will deal with the 
case of O bounded, it is straightforward to use our techniques to establish analogous 
well-posedness results for the problem posed on the whole space (i.e., H = K"). In 
fact, this problem is easier to approach as one has an exact representation of the 
associated heat kernel 



and this is the only significant alteration within the associated proofs. In the case 
of H = R", one would need to assume Lq, /q, and Vq decay suitably fast as jxj —>■ oo, 
for instance 


To(x) < C'|xr“ 


for some a > n, C > 0, and any x G H, in lieu of boundary conditions. Of 
course, for this particular application a bounded domain is certainly more natural, 
and we will typically have in mind the case of Neumann boundary conditions. 
Additionally, the cases of n = 1,2,3 are the most reasonable choices, but our 
results will actually be independent of dimension. Hence, we are able to derive sharp 
conditions on parameters that guarantee the local and global asymptotic stability 
of solutions regardless of the inherent spatial dimension. Additionally, though it 
is assumed throughout that A(x) is smooth, one may relax this assumption and 
arrive at similar conclusions regarding global-in-time solutions and their dynamical 
properties. Finally, another reasonable assumption to include within the model 
would be to take Dj = Dt, as infection of susceptible cells should not influence the 
rate of diffusion. However, such a condition will not be needed within the analysis, 
and hence we omit it. With the model well formulated, we state the main results 
of the paper. 


2.2. Existence of Solutions, Positivity, and a priori Bounds. Our first result 
is quite standard and serves to merely establish the existence of a suitable solution 
for small time. 
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Theorem 2.1 (Local Existence). Assume condition (Al) holds, then there is t* G 
[0, oo] such that (3CM) with (BCl)/(BC2) has a unique, positive solution (T, I, V) G 
C\iO,t*);C^n)) X Cl ((0,r);C2(ll)) X Cl {iO,t*);C\n)). 

We will omit the proof as the techniques are well-known in classical literature, 
and the result can be obtained by a straightforward application of the Contraction 
Mapping Principle. With a local-in-time solution in hand, we turn our attention 
to its properties, including positivity (assuming positive initial data), boundedness, 
and its extension globally in time. 

Theorem 2.2 (Positivity, Bounds, and Global Existence). Assume the initial data 
satisfy (Al). Then, for any t* > 0, there exist T, I, V satisfying (3CM) on (0,t*) x 
n, the initial conditions, and boundary condition (BC1)/(BC2). Furthermore, we 
have T(t,x),I(t,x),V{t,x) > 0 for all x € fl,t £ [0,t*], and there are C,Ct* > 0 
such that the solutions satisfy 

\\Tmo.<c 

and 

\\im^ + \\v{t)\\^<c,. 

for all t G [0, t*]. 

A few comments regarding the global existence of solutions are in order. Specih- 
cally, it should be noted that the idea of “mass transfer” among concentrations is a 
crucial aspect of this model. Consider, for instance, what would occur if the sign of 
the kTV term in the hrst equation were changed, so that T satisfied the evolution 
equation 

{dt - DtA)T = X{x) - htT -t kTV. 

Then, even in the simpler case of A = 0, it can be shown that this model mimics 
the behavior of the related nonlinear system 

{dtu — Au) = uv, {dtv — Av) = uv 

under a suitable transformation of the unknown quantities. Using tools similar to 
those established in [5, 32], it can be shown that solutions to this system blow-up in 
finite time. Hence, the exact transfer, rather than addition of the nonlinear terms 
is crucial to the existence theory, as is the continued positivity of solutions. In 
addition. Theorem 2.2 demonstrates further that the effect of a continuing supply 
of T-cells, given by A(x) > 0, does not lead to a blow-up in the system 

Now that we have established that solutions exist globally in time, we note that 
generally, the diffusion operator {dt — DA) has a smoothing effect on initial data, so 
we expect some gain in regularity in the concentrations. More specifically, assuming 
that the initial data is in Lfijl), we expect that solutions not only remain in 
but actually possess weak derivatives in this space as well. 


Theorem 2.3 (Regularity). Let t* > 0 be given. If T,I,V satisfy (3CM) on 
(0,U)xU andTo,Io,Vo G (H) with (BC1)/(BC2), t/ien VT(t, •), V/(t, •), VU(t, •) G 
(H) for all t G (0,t*). In fact, T{t, ■),l\t, ■),V{t, •) G iL"* (H) for all t G (0, t*) 
and m G N. 
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2.3. Asymptotic Behavior of Solutions. Finally, with the well-posedness of so¬ 
lutions understood, we turn our attention to their longtime asymptotic behavior. 
In particular, we study both the local and global stability properties of the unin¬ 
fected steady state, described below. Of course, any time-independent solution, i.e. 
a triple {T{x), I{x),V{x)), must satisfy the nonlinear elliptic system 


-DtAT = X{x) - ^J.TT - kTV 
-DiAI = kTV - ^il 
-DvAV = NfijI - fivV 


( 2 ) 


which may have many solutions for differing parameter values. In particular, we 
focus on the unique uninfected state Ec := (Too(x), 0,0), which represents the ulti¬ 
mate clearance of the virus, where T^o satisfies the linear boundary value problem 


— DtAToo = A(a;) — lirToa (3) 

for a; S n with boundary conditions given by (BC1)/(BC2). We first summarize 
some useful properties of the steady-state T-cell count. 


Theorem 2.4. The function T^o defined by (3) with (BC1)/(BC2) is C°° and 
satisfies 


0 < Too{x) < 




for every x G ft with T^o ^ 0. 


In order to determine the stability of the aforementioned equilibrium state, we 
are first led to study the system (2) linearized about namely 


-DtAT = \{x) - ^ltT - kT^{x)V 
—DjAI = kTao{x)V — fijl 
-DvAV = Nfiil - fivV 


( 4 ) 


with boundary conditions (BC1)/(BC2). Here, the first equation decouples from the 
last two and we need only study the PDFs describing I{x) and V{x). Accordingly, 
we define the second-order, linear self-adjoint operator 


£ = V ■ (HV) -f M{x) 


( 5 ) 


on the Hilbert space 

T-L = {(j) G X H^{ft) : (fi and 4>2 satisfy (BC1)/(BC2)} 


where 


D = diag{DI, Dy) and M{x) 


kToQ {T) 
N^i —^v 


( 6 ) 


We note that £ is a negative operator if —M (x) is a positive semi-definite matrix. 
Because fij > 0 this latter condition is implied by det(—M(a;)) > 0 for all x € H. 
After some algebra, this condition is exactly ||i?o||oo < 1 where 


Rq{x) 


NkToo{x) 


( 7 ) 


h-v 

The stability properties of Ec then depend crucially on the greatest real part of 
the spectrum of £, namely 


770 := sup{Re(77) : 77 G cr(A)}. 


(8) 
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Due to the compactness of the inverse of it is well-known that (j(C) consists only 
of eigenvalues and, when ordered, these form a decreasing sequence of real numbers 
diverging to — oo. Additionally, tjq can be uniquely identified using the associated 
Rayleigh quotient 





Hence, the sign of the parameter 770 will determine whether all eigenvalues are 
negative, ensuring that solutions which begin sufficiently close tend to the equilib¬ 
rium state, or whether a positive eigenvalue exists, thereby creating instabilities 
within the system. Prior to stating the stability and instability results for the vi¬ 
ral clearance state, we first show that under particular conditions it is the unique 
biologically-feasible solution of the steady state system ( 2 ). 

Theorem 2.5. If rjo <0 then the triple (roo(a;), 0,0) is the unique nonnegative 
solution of (2) with (BC1)/(BC2). 

Next, we characterize the local stability properties of the clearance state. In the 
spatially-homogeneous case, as in many other biological systems, it is known that 
a single parameter dictates the dynamical structure of solutions. With regards to 
(1), the quantity known as the basic reproductive ratio, defined by Rq := 
describes the entire long time behavior of the infection. In this case the spatially- 
independent T-cell count is Too = 7 ^, which means that Ro{x) as defined by (7) is 
exactly the basic reproductive ratio. Additionally, the associated principal eigen¬ 
value is 770 = i?o ~ 1- Though Rq for the spatially-heterogeneous system is now a 
function rather than a single value, we may still characterize the dynamics in terms 
of the single parameter rjQ. Hence, the following local stability theorem establishes 
a parameter regime that generalizes the spatially-homogeneous model, in which the 
longtime asymptotics are still shown to be determined only by the sign of 770 . 

Theorem 2.6. LetT,I,V satisfy (3CM) with (BCI)/(BC2). //770 < 0 then the 
viral clearance state Ec is locally asymptotically stable. Ifpo > 0 , then it is unstable. 

Finally, we study the global dynamics of this steady state. Our first global result 
shows that the T-cell population tends to Too{x) exponentially fast for any initial 
data and within any parameter regime. 


Theorem 2.7. Let T,I,V satisfy (3CM) with (BCI)/(BC2). Then, for every 
t > 0 


\\T{t) - TooWoo < \\To - Toollooe-'^^* 


In view of Theorems 2.5 and 2.6, we see that for 770 < 0, Ec is locally stable and 
the only nonnegative equilibrium. Hence, one might expect that it is also globally 
stable in this case. While we do not have a proof of this result assuming only the 
condition 770 < 0 , we can do so under a slightly more restrictive assumption. 


Theorem 2.8. LetT,I,V satisfy (3CM) with (BC1)/(BC2). //||i?o||oo < 1 then 
there are Co,a > 0 such that 


||/(t)||oo + ||R(t)||oo <Coe-“‘ 


for every t > 0. 




STEPHEN PANKAVICH AND CHRISTIAN PARKINSON 


Combining this result with Theorem 2.7 yields a sufficient condition under which 
the viral clearance state is a globally asymptotically stable equilibrium point of the 
system with exponential decay. Though it is beneficial to characterize the dynamics 
of solutions, it may be difficult to explicitly determine T^o for given parameter values 
and domain 17. Hence, our hnal result provides an additional sufficient condition 
that can be readily computed to describe the global asymptotics of solutions. 

Corollary 2.9. Let T, I, V satisfy (3CM) with (BC1)/(BC2). < 1, then 

the conclusion of Theorem 2.8 remains valid. 

We note that the longtime behavior results are independent of both the diffusion 
coefficients Dt, Dj, Dy > 0 and the dimension n = 1,2,3. In the next section we 
prove the main results of the paper. 


3. Proofs of Main Results. As previously mentioned, the proof of the first the¬ 
orem is a fairly standard application of the Contraction Mapping Principle, and 
hence we omit it. Prior to proving Theorem 2.2, we establish a few minor lemmas 
regarding the scalar diffusion equation that arise from the Maximum Principle. 

Lemma 3.1. Let D > 0 and uo{x) > 0 for all x G LI be given and assume that u 
satisfies 

{dt - DA)u = g{t,x), xel7,t>0'l . . 

u(0, x) = uo(a;). x € LI j 

If g{t, x) > 0 for all X G Ll,t > 0, then u{t, x) > 0 for all x G LI, t > 0. 

Additionally, if g{t, x) > 0 for all x G Ll,t > 0, then u(t, x) > 0 for all x G Ll,t > 0. 
That is, {dt — DA) preserves positivity on [0, oo) x 17. 

The following corollaries, which will be useful throughout the remainder of this 
section, are then immediate consequences of this result. 

Corollary 3.2. Ifu{t,x) satisfies (9) with D > 0, then 

||w(l)|loo < llwolloo + / \\ 9 ir)\\^dT, for all t G[0,oo). 

Jo 

Corollary 3.3. Assume that u(t,x) satisfies the differential inequality 
{dt — DA)u < g{t, x), x G Ll,t > 0 
u{x, 0) = uo(x), X G LI. 

Then u{t,x) satisfies the same inequality as in Corollary 3.2. That is, 

||m(<)IIoo ^ ll^olloo + / IlffMlloo^A for all t G[0,oo). 

Jo 

With these results in place, we may now prove the global existence and positivity 
theorem. 


Proof (Theorem 2.2). Since the initial data is positive and solutions are smooth 
(see Theorem 2.3), the functions T{t,x), I{t,x),V{t,x) must remain positive for 
all X G 17 and up to some time within the interval of existence [0,t*). Define the 
maximal time of positivity within this interval by 


Q = sup 


jiG [o,r) 


T{s,x),I{s,x),V{s,x) > 0 for all x G 17, s G 
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Estimating on this interval we find 

{dt - DtA)T + ^tT < A(x). 

Multiplying by an integrating factor yields 

(dt - DTA)[ef^^*T] < \(x)e>^^K 

By Corollary 3.3, we find 

||nOe'^"‘IL<l|7^olL+ / 

^0 

and upon integrating 

IIUi)lloo < \\To\L ^ (1 - ■ 

flT 

In particular, we have a uniform bound for ||7’(t)||^ on [0,(3], denoted by 

rM:=||rolL + ^. 

flT 

Next, we use the positivity of solutions and this uniform bound within the equa¬ 
tions for I and V. On the interval [0,(3], we have (dt — DjA) I < kT(t,x)V(t,x). 
Using Corollary 3.3, this implies 

||/(t)IL<||/olL+fc /‘||T(r)U(T)|Ldr 
Jo 

<II^oIL + ^T’m / ||r(r)|Ldr 

Jo 

<c(^i+j^\\v(T)\\^dTy (10) 

Using the positivity of I within the equation for U, we have 

(dt - DvA)V < N^il(t,x) 
for t G [0,(3]. Invoking Corollary 3.3 yields 

\\V(t)\\^<\\Vo\\^+Nf,i f\\I(r)\\^dr 

Jo 

<cy+jy\i(T)\\^dTy (n) 

Finally, define 

</>(i) = l|/Wlloo + lir(t)IL, tG[o,Q]. 

Adding (10) and (11), we see 

4’(t) < C ^1 -I- J (j)(T)dy , t e [0,(3]. 

By Gronwall’s Inequality, we can conclude that 

(t)(t)<Ce^\ tG[0,Q]. 

Thus, both ||/(t)||oo and ||U(t)l|j^ remain finite on [0,(3]. Therefore, solutions can 
be continued indefinitely, as long as they remain positive. 

Now, these bounds will imply the continuing positivity of solutions to time t = Q. 
Rearranging terms within the equation for I yields 

(dt - DiA)I + ^iil = kT(t, x)V(t, x)>0, tG [0, Q]. 
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Multiplying this equation by an integrating factor, we arrive at 
{dt-DiA) = ke^‘^T{t,x)V{t,x). 

Thus, we have a diffusion equation for e^^*I{t,x) with a nonnegative forcing term 
and positive initial conditions. By Lemma 3.1, we can conclude that a:) > 0 

for all t € [0, Q], and hence I{t, cc) > 0 for t € [0, Q]. Positivity of V follows in the 
same way since V satisfies 

{dt - DvA)V + fivV = x). 

and this becomes 

{dt - DyA) [e^^V] = Nnie>^''*I{t,x). 

Thus, by Lemma 3.1, e^''*V{x,t) > 0 for t G [0,(5] and V{t,x) > 0 for t € [0,(5]. 
Finally, for the T equation, we use the local-in-time bound on ||P(t)||oo to find 

{dt - DtA)T > a - mtT - kT{t,x) ||P(t)|U 

> A - ^lTT - Ce^^T 

Rearranging yields 

{dt - DtA)T + [^lT + Ce^*] T>\ 

and using the integrating factor ({t) = exp |/g (/ir + Ce^'^) (ir|, we find 

{dt-DTA){at)T}>xat) 

and C{0)T{0,x) = Tq{x). Thus, by Lemma 3.1, we have ({t)T{x,t) > 0, and hence 
T{t,x) > 0 for t G [0,(5]. Therefore, solutions remain strictly positive throughout 
the interval [0, Q] and this implies Q = t*, the maximal time of existence. Of course, 
since the solution remains bounded on this interval as well, we find t* = oo and the 
proof is complete. □ 

Next, we turn to the proofs of the regularity results. For brevity, we will prove 
the first statement only. However, the same technique can be applied to derive 
estimates inductively and prove the latter result without the introduction of new 
ideas. For additional background on the specifics of obtaining the second statement, 
we direct the reader to [22, 23] as the higher-order regularity stated here can be 
deduced by straightforwardly adapting the arguments in those previous works. 

Proof (Theorem 2.3). Beginning with the equation for T, we multiply by T and 
integrate over H to arrive at 

liWTml-DT [ TATdx= j \Tdx-^lT\\T{t)\\l-k f T^Vdx. 

t/ Q ^ Cl ^ Cl 

Integrating by parts on the left, enforcing boundary conditions, using the Cauchy- 
Schwarz inequality, and replacing V by it’s supremum, we see 

~\\T{t)\\l + DT\\VT{t)\\l 

< IIAII^ ||r(t)i|, - ||T(t)l|^ + fc||H(t)IL ||r(t)|i^ 

< \ (||A||^ + llT(t)ll^) - HT \\T{t)\\l + k ||H(t)|U l|T(t)||^ . 

Finally we arrive at 

I ||T(t)|l^ < a. (l + \\T{t)\\t) - 2 Dt\\ yT{t)\\l. 


( 12 ) 
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To generate a derivative estimate, we proceed similarly and take the gradient of the 
equation dotted with VT so that 

^dt I VT|^ - DtVT ■ VAT = VT ■ VA - fxr | VTp - kVT ■ (VVT + TVV). 
Integrating over spatial variables and using integration by parts yields 

~\\VTml + DrWATml 

= J VT-VXdx-fiT\\VT{t)\\l-kJ (v \VTf + TVT ■ VV^ dx 
< ^ (livr(t)ii^ + ||va|| 2) - MT l|VT(t)||2 + k (lini)IL w^mwl + 

+ ^\\m\\^[\\VTit)\\l + \\VV{t)\\l)), 

whence 

^ ||vr(t)||^ < a. (l + ||vr(t)||^ + ||y(t)||^) - 2 Dt\\AT ml ■ (13) 

We deal with / analogously. Multiplying this evolution equation by I and integrat¬ 
ing gives 

^ wnt) wl+Di 11 v/(t) 1 12 = ^ ^ Tivdx - \\m 11 ^ 

<lini)IL / TIdx-fij\\m\\l 

Jn 

<h\\vit)\\^{\\Tml + \\m\\l)-t^i\\iml 

From this we find 

I wmwi < o. (iir(f)ii^ + ii/(t)ii^) - 2 Dj wvmwi. (m) 

Next, taking the gradient of this equation and then dotting with VI, we arrive at 

^dt \VI\^ - DiVI ■ VAI = kVI ■ {VVT + TVV) - \VI\^. 

Integrating over ft yields 

~\\VI{ml + Di\\AI{t)\\l+ ^^J\\VIit)\\l 

= k(^J VVI-VTdx + J TVI-VVdx'^ 

<k(\\V{t)\\^[ VI-VTdx+\\T{t)\\^ [ Vl-VVdx) 

\ JQ Jn J 

< c*. (^1 (||v/(t)||^ + \\vT{t)\\l) + i (||v/(t)||^ + livy(t)ii^) 

and this produces the inequality 

I ||V/(t)||^ < Ct. {\\VT{t)\\l + wvmwl + \\VV{t)\\l) - 2Di\\AIit)\\l. (15) 
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Within the V equation we multiply by V and integrate to find 

ii^wii2+ w^yml = ^ ivdx - f,v\\v{t)\\i 

<^{\\iml + \\vml)-dv\\vml 
<Ct^ (ii/wii2 + iini)ii2)- 

This yields 

I \\vml < ct^ {wmwl + ii^wii") - ‘ 2 Dv wvvml ■ (le) 

The analogous derivative estimate is then 

W^vml + Dv WAV ml = VI ■ VVdx - W^vml 

< ^ (liv/(t)||^ + livy(t)ii^) - ^^v W^vml 

<Ct^ (||V/(t)||^ + ||Vl/(t)||^), 

from which it follows that 

I w^vml < m (iiv/(t)ii^ + iivF(t)ii^) - 2Dv wAvml ■ ( 17 ) 

Finally, let Z?min = inii[i{DT, Di, Dy} and for t S define 

M{t) = [wTml + wnml + wv{t)\\i) 

+ ||VT(t )||2 + ||V/(t)||^ + WVVml ) 

=: 00(t) +t4>i{t). 

By adding equations (12),(14) and (16), we see 

0o(O < Cy (1 + 0o(i)) - 2D^im(t) 

and by adding equations (13), (15) and (17), we see 

Mt) < a.(l + 01 (t)) - ||AT(t)||^ + ||A/(t)||^ + WAVml ) ■ 

Then, letting 02(t) = ||AT (<)||2 + ||A/(t )||2 + ||A14(t)||2 , we arrive at 
M\t) = 0'(t) + “t” -^min 

C-f^* (1 -j- 0Q(t)) 2Z?jnin0l(O ^min^ (^Ct.(l + 0l(<))- 202(t)) 

< Ct.(l + 00(t) + ^min t0l(t)) - -71min0l(t) 2Zlniint02 (t) ■ 

Noting that 0i and 02 are nonnegative, we find 

M'{t)<Cy{l + M{t)). 

From this, an application of Gronwall’s inequality produces the bound 
M{t) < Cf (1 + M(0)e*) < Ct* (1 + M(0)). 

By assumption 

M(0) = ||ro||^ + ||/o||^ + 11^0112 
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is finite and thus M (t) remains finite on the interval. Since (j)o is nonnegative, this 
implies the bound 

for t G (0,t*]. As each of the quantities ||Vr(t )||2 , ||V/(t )||2 , ||Vy(t )||2 is bounded 
by we find 

vr(t,-), vv{t,-)GLAn) 

for any t G (0,t*]. To prove the latter statement of the theorem, higher-order 
estimates are obtained and combined using the same techniques. The result then 
follows by induction (see [23] for further details). □ 


With the regularity proof complete, we finally study the longtime dynamical 
behavior of (3CM), and in particular, the viral clearance steady state -Ec- We first 
prove some useful properties of Too (x). 


Proof (Theorem 2.4)- To prove the nonnegativity of Too, we merely note that the 
Green’s function associated to the operator (—-b fix) is positive [14], and 
hence the conclusion follows due to the nonnegativity of A(a;) on fl. Additionally, 
Too inherits the regularity of A and Too A 0 because G, A ^ 0. 

To prove the upper bound, we let G{x) satisfy 

4 - iit)G = 5{x) 


for X G Tl. Then, integrating over fl, we find 

f , 1 Dt f dG 

/ G{x) dx = - 1 -/ dS. 

Jn MT Jan on 

If G satisfies (BC2), we merely enforce the boundary conditions to conclude 

f G{x) dx = 

Jn MT 

If G satisfies (BCl), we note that G = 0 on dil and G{x) > 0 for a: G fl, therefore 
^ < 0 on dfl and thus 

on — 

/ G(x) dx < —. 

Jn Mt 

Again, since G is positive, this further shows UGjji < ^ in either case. Finally, 
since Too satisfies the inhomogeneous equation it must be given by Too = G * A and 
a simple convolution estimate yields 


[Toolloo = ||G*Alloo < 


.||Glli = 




□ 


Next, we prove the uniqueness of this steady state when rjo < 0. 

Proof (Theorem 2.5). Let {T{x), I{x),V{x)) ^ (Tao(x),0,0) be another nonnega¬ 
tive solution of (2) satisfying (BC1)/(BC2). In particular, we must have / ^ 0 and 
V ^ 0 since either condition implies the other and T = Too. The second and third 
equations in the system, namely 

-T>jAI = kfV - nJ 
-DvAV = NfiJ - fiyV 
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can be rewritten as 


or 


DiAi + kTooV - fill = k{Too - f)V 
DvAV + N^ii - nvV = 0 


/ ■ 


■ k(Too - f)V ■ 

_ V _ 


0 


(18) 

Next, let u{x) = T(x) — Too (a;) where T^o is defined as the unique solution of the 
linear equation (3) with (BC1)/(BC2). Then, u satisfies 

[—DtA + /It] u = —kTV. 

Since T{x), V(x) > 0 for all a; G O, it follows that u{x) < 0 for all a; G O. Therefore, 
equation (18) can be rewritten as 

C 


and integrating over 11 , we find 


I ■ 


—kuV 

V 


0 


Taking the dot product of this equation with 


I 

V 


f 

i{x) 

■ c 

i{x) 

dx = k I 

Jn 




Jn 


However, since rjQ < 0 and I,V ^ 0, we find 


I{x) 


■C 


_ V{x) \ [ V{x) _ 


i{x) 


dx < rjo 


2 + 


< 0 . 


This inequality contradicts the previous one, and hence the original assumption of 
an additional solution cannot hold. □ 

Next, we prove the local stability and instability results. Since a generalization of 
the celebrated next generation method adapted to spatially-heterogeneous problems 
was recently developed in [33], we will rely on many of the results therein for the 
proof. 

Proof (Theorem 2.6). To prove the first conclusion, we assume 770 < 0 and put the 
system (4) into the form of [33, Equation (3.9)]. The uninfected compartment - 
T{t, x) - and infected compartments - I(t, x) and H(t, x) - can be separated into the 
scalar function us and vector function u/, respectively. The equation for the latter 
can be represented as 


dtui = V ■ {DVui) + F{x)ui — V{x)uj 


(19) 


where 


and 


ui{t,x) = 


I{t,x) 

V{t,x) 


D = diag{Di,Dv), 


F{x) = 


kTaa{x) 

0 


V{x) = 


Ti 


0 

Tv 


Using the properties of Too guaranteed by Theorem 2.4, we note that the system in 
this form satisfies assumptions (A1)-(A6) of [33]. Then, [33, Theorem 3.1] directly 
implies the conclusion under the assumption 779 < 0, since this eigenvalue is exactly 
the spectral bound for the operator on the right side of (19). Therefore, is a 
locally asymptotically stable equilibrium for (3CM). 
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To prove the second conclusion, we will show that the eigenvalue problem 

V • {DVcj}) + = r](j) (20) 

with D and M defined by (6) possesses a positive eigenvalue with corresponding 
positive eigenvector under the assumption r/g > 0. To this end, we will again utilize a 
result from [33]. From (6), we see that M(x) is cooperative for all a; G 12 and M{xg) 
is irreducible for some xq G since, by Theorem 2.4, Too{x) > 0 for every x G 12 
and Too ^ 0. Hence, by [33, Theorem 2.2], rjo as defined by (8) is an algebraically 
simple eigenvalue of (20) with a strongly positive eigenvector satisfying 

Re( 77 ) < 770 for all 77 G cr^V • {DV) + \ { 770 }. 

Therefore, 770 is a positive eigenvalue with positive eigenfunction, and the instability 
of Ec follows immediately. Finally, we note that that these results are valid for either 
boundary condition (BCl) or (BC2) as mentioned in [33, Remarks 2.2 and 3.1]. □ 

With the local stability proof complete, we turn to the proofs of the global 
asymptotic behavior of the system. 

Proof (Theorem 2.7). Similar to the proof of Theorem 2.5, we define u{t,x) = 
T{t,x) — Too{x) so that u satisfies 

dtu — DtAu = —^tu — kTV. 

Using the integrating factor and recalling the positivity of T{t,x) and V{t,x) 
guaranteed by Theorem 2.2, this becomes 

{dt - £>TA)[e'"^*u] = < 0. 

Invoking Corollary 3.3 with g = 0 yields 

l|e'"^‘M(t)l|oo < ||m(0)11oo 

which is equivalent to 

||T(2) - Toolloo < llTo - Toollooe-'^^* 

for every t > 0, and the proof is complete. □ 


Proof (Theorem 2.8). We begin by writing the system in non-dimensionalized form. 
Although one can arrive at the same results without doing so, using the dimension¬ 
less system will simplify the proof by reducing the number of parameters which 
appear and the complexity of related expressions. Only within this proof will the 
non-dimensionalized version of the system be utilized. 

First, define the dimensionless concentrations by 


T*{t,x) 


T{t,x) 

Tc 


r{t,x) 


Ic 


V*{t,x) 


V{t,x) 

Vc 


where Tc^Ic, and 14 are constants to be determined. Additionally, we scale the 
spatial and time dimensions by letting 


t* = 


t 

4’ 


X 

Xc 


Here the scaling constants will be determined so as to minimize the dimension of 
the resulting parameter space. Of course, differing spatial parameters (i.e., one in 
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each spatial component) could be introduced, but we assume a uniform directional 
scaling. Substituting these expressions within (3CM), we find 


dT* 

~ 

dr 

W 

dV* 


DricAT* = ^x{x) - ^irtcT* - ktcV,T*V* 


dt* 


xi 

_ 

Dytc 


Tc 

■AI* = k- 




AV* -fivtcV*. 

y r. 


( 21 ) 


We fix the time and spatial scales using the T-cell diffusion and decay rates by 
choosing tc = and Xc = \J Next, we choose the scaling for the dependent 
variables so as to eliminate parameters in each equation. In particular, this is 
effective when 


T,= 




Ir = 


fJ-VfJ-T 


Vc = 


fJ-T 


kN’ " kNm' " k 
Dropping the starred notation, this finally yields the non-dimensionalized system 

dT 


dt 

dt 


- AT = q{x) -T -TV 


- = ai{TV - I) 

- P 2 AV = a2{I -V) 


( 22 ) 


where 


and 


Oil = 


tL 

flT ’ 


k-V o Di 

Oi2 = -, PI = 7 ^, 

flT 


R 


q{x) := 


kN 


X{x). 


kTkv 

We note that the system now possesses only five parameters, and as we will see, 
the a terms control its rates of convergence to equilibrium, while the /3 coefficients 
describe the rates of relative diffusion. Finally, we rescale the steady state Toc{x), 
defined by (3), using the same values of Xc and Tc so that, again dropping the 
starred notation, this function satisfies the linear, stationary PDE 


- AToo = q{x) - To, 


(23) 


Due to the introduction of the scaling parameter Tc, our initial assumption of 
lll?o||oo < 1 then implies the newly scaled equilibrium T-cell count satisfies || Tool! 00 < 
1. Recall that all original parameters are assumed to be positive, and hence all new 
parameters retain this property. To prove the result, then, we will show that the 
conclusion holds for solutions of (22) as the exponential decay of the quantities / 
and V in the scaled space-time variables will imply the result of the theorem in the 
unsealed variables with a change in the constants Cq and a. 

With the dimensionless equations in place, estimates for T are immediate. From 
Theorem 2.7 we rescale the resulting inequality (or equivalently, perform the same 
steps within the proof on (22) and (23)) to find 

||r(t)||oo < llToolU + l|r(0) - Toollooe-* =: Pit). 
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As ||Too||oo < 1, let us define the quantities M = ^(1 + ||T'oo||oo) and M = i(l + M) 
so that ||Too||oo < M < M < 1. Since P{t) is decreasing, we further choose 

ro-To,ii< 


To = max < 0, In 


M - \\T^ 


so that P(t) < M for all t > tq. 

Next, let a* = ^ inin{l, ai, 02 } and choose a € (0, a*) small enough such that 

Q;ia2 


b:=M - 


(ai — a){a 2 — a) 


M 


satisfies & > 0. This choice of a will allow us to take the constant in the exponential 
decay sufficiently large so that T is uniformly bounded and the decay is preserved 
in time. In particular, let Ci > 0 satisfy 


^ . .fsup,g[ 0 ,.„] ||I"(s)||oo SUp,g[o,ro] IIAs)lloo] 

■ I T(o)iu ■ iino®: r 

d 

^ ^ 1 / Cl2||J(0)||oo \ 

^-b \ ^ {a2-a)\\V{0)U)' 

With these constants in place, notice that for every t G [0,to], we have 
imoiloo < sup |iy(s)lloo < Ci||y(0)||ooe-“"“ < Ci||F(0)|Ue-“‘. 

sG[0,ro] 


(24) 

(25) 


Thus, the definition of Ci implies that V satisfies the decay estimate on the bounded 
interval [0,to], and the same inequality holds for /. It remains to prove these 
estimates for t > tq. 

Now, because solutions are sufficiently regular we may invoke Theorem 2.3 and 
the Sobolev Embedding Theorem to deduce that the mapping t —>• ||E(t)||oo is 
continuous. Hence, we will utilize a continuity argument involving this norm in 
order to complete the proof. Let 


Ti = sup{t > 0 ; ||H(s)||oo < ^ 1111 /( 0)11006 for every s G [0,t]} 
and notice that ti > tq. Then, for t G [ro,Ti], we have both 

||H(t)||oo <Ci||I/(0)||ooe-“‘ and ||T(t)||oo < M. (26) 

Deriving estimates on I{t,x) as for its dimensionalized version, we see 

- /3iA^ [e“i‘/] = aie“i*rH. 

Integrating and using both Corollary 3.2 and (26) implies 

||/(i)||oo < ||/(0)||ooe-“^‘ + aiMe-“i* / Ci||I/(0)||ooe(“^-“)* ds 

< |l/(0)||ooe-“^‘ + ^^||I/(0)||ooCiM (e-“‘ - 6-“^*) 

Qi — a 

< C2e-“‘ 

where 

C 2 = Il/(0)||oo + ^^||H(0)||ooCiM. 
ai — a 

Doing the same for V and using the newly-derived estimate on ||I(i)||oo, we find 

- /32A^ [e“=V] = a 2 e“^‘J 
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and thus 

^0 

< imO)||ooe-“^‘ + ^^C'2(e-“‘-e-“=‘) 

a2 — a ^ 

< f||nO)||oo + ^^C2)e-“‘. 

V a2-a J 

Thus, the exponential decay of ||l^(t)||oo continues on this time interval and we 
merely require the constant within this inequality to be strictly dominated by 
Ci||y(0)||oo to complete the argument. After a brief calculation, we see that (24) 
implies 

iino)iioo + -^^C2 < MCiWvmu 

a2 — a 

and combining this with the estimate above, we find 

||f"(i)||oo<MCi||F(0)||ooe-“‘ 

for any t G [To,ri]. Since M < 1, we see that ti cannot be finite, as this would 
contradict its definition as the supremum, and thus ti = c». Therefore, using (25) 
the estimate 

|ini)lloo<Ci||F(0)||ooe-“‘ 

holds for all t > 0. Since we have ||/(t)||oo < C 2 e~°‘* on the same interval, the 
exponential decay of ||/(t)||oo follows with the same rate. Taking Co = maxjCi, C 2 } 
yields a uniform decay estimate on both non-dimensionalized quantities, namely 

\\im^ + \\vmoo<c,e-^'. 

Finally, upon rescaling the non-dimensionalized system, the result holds for the 
original concentrations I and V with a change in the constants Co and a. 

□ 


We end this section with the proof of the associated corollary. 


Proof (Corollary 2.9). Assuming the condition on the parameters holds, we note 
that ||Too||oo < by Theorem 2.4. Hence, < 1 implies 


||d?o||oo 


A^fcllT’oolloo 


< 1, 


fJ-V 

and the assumptions of Theorem 2.8 are satisfied, which implies the result. 


□ 


4. Proofs of Lemmas. Finally, to complete the paper, we include the proofs of 
lemmas from the previous sections. The first few results are fairly straightforward 
applications of the ideas inherent within the Maximum Principle, but we include 
them for completeness. 

Proof (Lemma 3.1). While both results are classical, we will prove the former asser¬ 
tion and direct the reader to [6] for details regarding the latter. Define the positive 
and negative parts of u by 

(t, x) = max{0, u(t,x)} 

U-(t,x) = — min{0, u(t, x)} 
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and notice that u = u+ — u-. With this in mind, we multiply (9) by U- and 
integrate in both time and space. Then the left side is given by 


L ■= 


U-dtU dxdt — D 


0 Jn 


u-Au dxdt ■■= I + II, 


0 Jn 


where s > 0. Now, define = {{t,x) G (0,s] x : u{t,x) < 0}, then it_ = 0 
outside of 0 “ so we see 

I = u-dtudxdt. 

J Q,s 

However, on this set u = —u_ so 

I = — j U-dtU- dxdt = — 7 : [ dt (w?-) dxdt. 

However, since u_ = 0 outside of Hj, integrating over Hj is the same as integrating 
over (0, s] x H. Thus 

/ = —i f [ dt {ut.) dtdx =—\- f [u-{s,x)‘^ — u-{0,x)'^] dx. 

^ Jn Jo 2 

Finally, uq > 0 implies that U-{0,x) = 0 so we find I = —^ ||u_(s )||2 < 0. 

We use a similar set of steps for II: 


II = -D 

= D I 

= -D 


U-Au dxdt 

lo JQ 

U-Au- dxdt 


'K 


Vu_ • Vit- dxdt + D u- 


du- 


dn 


dt. 


dn 


Enforcing the homogeneous boundary condition, either (BCl) or (BC2), yields 


II = -D 


n \Vu-{t,x)f dxdt < 0. 
! 


Thus we have 

L = / + //< 0 . 

Next, considering the right side under the same operations, we have 


R := 


:= / g(t,x)u-{t,x) dxdt. 

Jo Jn 


But both g and U- are nonnegative so i? > 0. 

Because L = R with L < 0 and i? > 0, we find that both must be zero. Addi¬ 
tionally, L = 0 forces I = II = 0, and we conclude 


I«-W||2 


= 0 . 


This is only possible if u_(s, x) = 0. However, s was an arbitrary element of (0, 00 ) 
and so U-{t, x) = 0 for all x G fl, t G (0, 00 ). But if u- = 0, then u = it+ > 0 which 
completes the proof. □ 


Proof (Corollaries 3.2 and 3.3). Define the function 
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Then we notice that 

{dt-DA)v = ||g(t)|L - {dt-DA)u = HffWIL “ 2 ;) > 0. 

Also 

n(0,a;) = ||uo|L “ ^(0:a;) = HmoIL “ > 0 . 

Thus, by Lemma 3.1, v{t,x) > 0 for a: € n, t > 0, and 

u(t,x) <\\uo\\^+ [ \\g{T)\\^dT, 

Jo 

for all X € r2, i > 0. Taking the supremum over x € ft yields 

IImWIL < ll^olloo + / \\9{T)\\^dT, 

Jo 

for all i > 0 which completes the proof. The same method then applies to prove 
Corollary 3.3. □ 
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